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Abstract 

Two-particle correlation functions of negative hadrons over wide phase space, and trans- 
verse mass spectra of negative hadrons and deuterons near mid-rapidity have been mea- 
sured in central Pb+Pb collisions at 158 GeV per nucleon by the NA49 experiment at 
the CERN SPS. A novel Coulomb correction procedure for the negative two-particle cor- 
relations is employed making use of the measured oppositely charged particle correlation. 
Within an expanding source scenario these results are used to extract the dynamic char- 
acteristics of the hadronic source, resolving the ambiguities between the temperature and 
transverse expansion velocity of the source, that are unavoidable when single and two 
particle spectra are analysed separately. The source shape, the total duration of the 
source expansion, the duration of particle emission, the freeze-out temperature and the 
longitudinal and transverse expansion velocities are deduced. 



1. Introduction 

Lattice QCD calculations |I] predict that a phase transition occurs between hadronic 
and deconfined, chirally symmetric partonic matter at energy densities of about 2 GeV/fm 3 , 
corresponding to a critical temperature T c between 150 and 200 MeV in baryon-free mat- 
ter. Studies of transverse energy distributions in central Pb+Pb collisions at the CERN- 
SPS with 158 GeV per nucleon 208 Pb beams (y^s « 17 GeV per nucleon pair) show that 
energy densities of 2-3 GeV/fm 3 are created in the initial interaction volume @. It is 
therefore conceivable that such a transition is induced in ultrarelativistic nuclear colli- 
sions. The partonic phase created in the initial stage of the collision hadronizes as the 
system expands and cools. The expansion dynamics of the final hadronic phase are ex- 
pected to differ considerably with the presence or absence of a prior plasma phase, and it 
is suggested that the observed final state hadrons reflect the space-time evolution of the 
system and therefore provide clues as to the state of the system prior to hadronization. 

The space-time evolution of a pion emitting source can be probed using the momentum- 
space correlation of identical pions. Bose-Einstein correlations have been first used to 
study the properties of a hadron emitting system created in hadronic collisions under the 
assumption of a static source JJ. The expansion of a source that is not static but exhibits 
a hydrodynamical scaling evolution has been discussed by Shuryak and Bjorken 
Within the framework of such models it has been further shown that the dependence of 
the correlation function on the pion pair kinematic variables is an important diagnostic 
tool for the study of an expanding system |j|-|8j . 

Experiments NA35 |9| and NA44 |nj have applied this technique to the investigation of 



central collisions of 32 S nuclei with heavy nuclear targets at the CERN SPS, resulting in the 
picture of a cylindrically symmetric source expanding predominantly in the longitudinal 
direction, with some indication of concurrent transverse expansion. NA35 has observed 



that the incoming energy is not completely stopped for these systems [EC , consistent with 
the predominance of longitudinal expansion. The longitudinal evolution was found to be 
consistent with Makhlin and Sinyukov's analysis |J of "scaling" expansion, which is based 
on Bjorken's earlier study || of an idealised hydrodynamic longitudinal expansion model. 
Longitudinal collective expansion is not unique to relativistic nucleus-nucleus collisions. 
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However the transverse expansion dynamics, seen through single particle and two particle 
observables, is expected to be sensitive to the existence of a first order hadronization 



transition 0, as well as to the occurrence of a near isentropic hadronic expansion |12 
which would preserve information about the initial phase, and the existence of a primordial 
partonic transverse expansion velocity fll3yi4| . 

Transverse collective expansion is expected to be more prominent in the heavy Pb+Pb 
collision system than in collisions of lighter projectiles, owing to the smaller surface-to- 
volume ratio, the larger "fireball" volume and the longer duration of expansion. Indication 
that this is the case is provided by the single particle transverse mass spectra of kaons, pro- 
tons and As produced in central Pb+Pb collisions at the CERN SPS Hl5,16|, which exhibit 
unusually high inverse slope parameters ("temperature"), in excess of 250 MeV, ruling 



out a simple thermal hadronic gas model |L7|. These distributions can be reproduced by 



invoking a radial velocity field in transverse direction, though the single particle spectra 
alone cannot distinguish between an elevated temperature and a boost due to transverse 
collective expansion [|T^,|T^] and therefore cannot determine uniquely the temperature and 
transverse velocity. 

It will be shown that the combination of single particle spectra with two-particle corre- 
lation results can lift this essential ambiguity and determine both the expansion velocity 
and temperature at freezeout. The term freezeout refers here to the final decoupling of 
hadrons from the strong interaction. Note that the term "freezeout" is often also em- 
ployed for hadronic decoupling from inelastic interactions at which point the population 



ratios of hadronic species become stationary |[9| . This stage obviously precedes the one 
under study here. 

The analysis of this paper is based upon measurements by the NA49 experiment of 
two-particle correlation functions of negative hadrons over a wide phase-space region and 
on single particle transverse mass spectra of negative hadrons and deuterons near and 
below mid-rapidity, for central Pb+Pb collisions at 158 GeV per nucleon. The like-charge 
pion correlation functions are corrected for Coulomb effects using a novel technique based 
upon the parametrization of unlike-charge pion correlation functions measured in the same 



experiment [2(J . The corrected like-charge correlation functions are fitted with the Yano- 



Koonin-Podgoretskii parametrization [21J as formulated by Heinz and the extracted 
source parameters are studied as function of rapidity and transverse momentum of the 
particle pair. In the framework of an expanding hadronic source model, such an analysis 
of correlation functions yields information about the source size and the duration of the 
freezeout process, as well as the ratio /3j_/T near mid-rapidity, where f3± is the transverse 
velocity and T is the freezeout temperature. The ambiguity between temperature and 
expansion velocity implicit in the ratio can then be lifted by considering single particle 
transverse mass spectra. A consistent picture of collective velocity fields governing lon- 
gitudinal and transverse expansion emerges, with an estimated freezeout temperature at 
mid-rapidity of 120 MeV and a transverse expansion velocity of f3±_ = 0.55. 

The paper is organized as follows: section || describes the experimental procedure and 
data analysis, section |^ gives a discussion of the experimental determination of correlation 
functions, section |4] describes the determination of the source parameters through fitting 
with the Yano-Koonin-Podgoretskii parametrization of the correlation functions, section 
D discusses the extracted parameters in terms of source size, average freezeout time and 
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Figure 1. Schematic diagram of the NA49 experiment. The beam entering from the left 
interacts with the target T, producing charged particles most of which are analyzed by 
two superconducting magnets and their trajectories are recorded by four Time Projection 
Chambers (TPCs) with an active volume of 50 m 3 . Particle identification is done by time 
of flight and dE/dx. A set of calorimeters cover the forward rapidity region and are used 
for triggering purposes. 



duration, freezeout temperature and expansion velocity, section § introduces transverse 
mass spectra to determine uniquely the freezeout temperature and expansion velocity 
within the framework of a hydrodynamical model, and conclusions are drawn in section [7[ 

2. Experimental Procedure and Data Analysis 



The NA49 experiment [ IB| , shown schematically in Fig. ^ is a large acceptance hadron 
spectrometer at the CERN SPS. It comprises four large Time Projection Chambers 
(TPCs), four Time-of-Flight walls (TOF), and a forward (Veto) calorimeter for trigger- 
ing. The target (T) was a foil of natural lead with thickness 224 mg/cm 2 (1% interaction 
probability for the 208 Pb beam ), mounted 80 cm upstream of the first TPC. 

Two of the TPCs (VTPC1 and VTPC2) are placed inside dipole magnets of 9 Tm 
total bending power and achieve a momentum resolution 5p/p 2 = 10 -4 GeV -1 . The 
TPCs downstream of the magnets (MTPCs), extend the acceptance to high momentum 
and provide particle identification via ionization measurement with a resolution of 4.5% in 
dE/dx. The rapidity coverage of the TPCs is 2 < y n < 5.5, and the transverse momentum 
range selected for this analysis is 0.05 < px < 2 GeV/c, with almost complete azimuthal 
coverage. The Time-of-Flight walls cover a rapidity of 2 < y^ < 3.5 with an average time 
resolution of about 60 ps and a flight path of about 14 m. 

The particle identification using the measurement of the specific ionization in the TPC 
gas in the relativistic rise regime is complemented by the time-of-flight measurement at 
the lower momenta. In the present work we utilize particle identification only for the 
deuteron transverse mass spectra. Two-particle correlations and pion transverse mass 
spectra are based upon measurement of negative hadrons, which are dominantly pions 
(see below). 

A prompt trigger selects events from the lowest 3% of the forward energy distribution 
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detected in the Veto calorimeter, located 27 m downstream of the target and covering the 
beam fragmentation region (laboratory rapidity y > 5.5). This anticoincidence trigger, 
which vetoes events with a large number of projectile spectator nucleons, serves as a 
"central collision" trigger selecting impact parameters less than 3 fm [Q. About 1,200 
charged particles are recorded in the NA49 TPCs for each central Pb+Pb collision. 

Negative hadron two-particle correlation functions are studied using 80,000 events an- 
alyzed independently in the second Vertex TPC and in the Main TPCs. Negative hadron 
and deuteron transverse mass spectra are analyzed using 300,000 events in which charged 
particles are tracked from the MTPCs to the TOF walls. The detailed techniques em- 
ployed for cluster finding, tracking, efficiency calibrations, non-vertex track rejection, and 
hadron identification |23| will be described in a forthcoming publication. 



Since the two-particle correlation function carries the information of the source size in 
the small relative pair momentum range, and this region becomes narrower in momentum 
space the larger the source size is, the relevant signal for central Pb+Pb collisions is 
expected to occur at momentum differences less than 50 MeV/c in the pair centre of 
mass frame. Therefore the required resolution for the measurement of the momentum 
difference of two tracks is less than about 10 MeV/c. The momentum resolution was 



studied using a Monte Carlo simulation chain based on the event generator VENUS |24| 
and the detector response simulation code GEANT |p5| . The resolution in all relative 
momentum projections was shown to be (5 ±1.5) MeV/c throughout the acceptance used 
in this study, 2.9 < Y n7T < 5.5 and pair mean transverse momentum K± < 0.6 GeV/c 
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3. Discussion of Correlation Functions 

Experimentally, the normalized two-particle correlation function is obtained from the 
ratio 

where A(Q) is the measured two particle probability distribution as a function of relative 
four momentum Q = pi — p 2 , and B{Q) is the (uncorrelated) background distribution, 
calculated in the same way but using two particles taken from different events. 

The dynamical expansion of the source in central collisions is reflected in the depen- 
dence of the correlation function on all five non-trivial momentum components of the 
pair kinematics, in which the azimuthal angle with respect to the beam axis is averaged 
out. For extraction of the source parameters from the correlation functions we employ 
the Gaussian approximation in a representation given by Heinz that is appropriate to 



explicitly reveal longitudinal collective motion [p,22 : 



C 2 (Q±,Q\\,Qo;Y nn ,K ± ) = l + Xexpl-QlRl-^iQw-PYKQofRl- 

Jyk(Qo-(3ykQ\\) 2 R 2 ]. (2) 

This expression is an explicit function of the energy difference Qo = Ex — E 2 , trans- 
verse momentum difference Q± = ((Ap x ) 2 + (Ap^) 2 ) 1 / 2 , longitudinal momentum differ- 
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ence Qu = Ap Z} and correlation intensity parameter A, and an implicit function of the 
pair rapidity 

y^ = Ih e J + e J +Pz1+Pz2 ). (3) 

2 K E 1 + E 2 -p gl -pa J 1 ; 



and the mean transverse momentum K± = \\J(p x i + Vx2) 2 + (p y i +P y 2) 2 - Pyk is the 
"Yano-Koonin velocity" which describes the source's longitudinal collective motion in 
each interval of Y n7T . A non-zero (5yk results in mixing the longitudinal momentum and 
energy differences, Q\\ and Qo, in determining the longitudinal source parameter Ru and 
the time duration parameter R . The parameters A, R±, R\\, R and Pyk are extracted 
from a fit to the correlation function data. The parameters R±, Ru and Rq describe 
the apparent source extent in transverse, longitudinal and temporal directions. They 
still depend on pair rapidity Y n7T and pair transverse momentum K± . This dependence 
will be further analysed in sections 4 and 5, resulting in the extraction of the actual 
geometric space-time parameters of the source within an explicit model with longitudinal 
and transverse expansion ||. 

Due to its large acceptance, NA49 was able to determine the dependence of expres- 
sion (pj) on all five kinematic quantities over a wide region of Y n7T and K± , giving a 
complete picture of the source expansion dynamics. 

To investigate the effect of the finite resolution and other experimental effects on the 
correlation functions, a simulation of the Bose-Einstein correlations |2B| was applied to the 



(uncorrelated) VENUS events, which were then processed through the detector simulation 
and the same analysis chain that was applied to the data. In this way, the response of the 
device was estimated to generate a lowering of the source parameters i?o, R± and R\\ by 
an amount between 2 and 8%. No correction for this effect was applied to the parameters 
reported below. 

The correlation functions were calculated using all negative hadrons rather than identi- 
fied 7T~ pairs. The Monte Carlo simulation chain showed that the h~h~ pair signal consists 
of about 55% genuine 7r~7r _ pairs from the primary event vertex, with this fraction uni- 
form over the measured phase-space. The remaining background pairs arise principally 
from three sources: (i) K~, electron and antiproton contamination from the primary in- 
teraction; (ii) 7T _ from K and A decays in flight that were erroneously flagged as primary 
vertex tracks, and (iii) secondary tracks generated in the detector material that could 
not be removed by primary vertex cuts. However, simulation of the correlation functions 
comparing the h~ pair signal to the genuine ir~ pair signal showed that the contamination 
is "well behaved" , in that it contributes uniformly to all two-track relative momenta in all 
projections of the correlation function. As a result, the background contamination prin- 
cipally affects the absolute magnitude of the correlation signal but not the shape of the 
correlation functions. In an analysis employing equation (0), the Monte Carlo simulation 
thus showed that only the A parameter, which quantifies the overall correlation strength 
falls from its ideal value of unity down to A ~ 0.5 due to contamination in the true ir~ 
pair signal. The source parameters Rq, R± and R\\, which depend upon the shape of the 
correlation function, deviate at the level of only 2 to 6%. No correction for this effect was 
applied to the parameters reported below. Due to the large acceptance, high multiplicity 
and large number of analysed events, the statistical errors on the correlation functions are 
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negligible. Taking into account the systematical error of the fit procedures, we estimate 
the systematic errors of space-time parameters extracted from the correlation function to 
be ±15%. 

In order to extract the correlation function for identical charged pions which is due to 
quantum statistical effects, and which is sensitive to the characteristics of the emitting 
source, it is necessary first to correct the observed correlation for the effect of Coulomb 
repulsion. The distribution in relative momentum of two charged pions is strongly affected 
by Coulomb repulsion as their momentum difference, Q, approaches zero. The correction 
for the Coulomb repulsion is done by applying a weight to each pion pair in the background 
distribution B(Q) in equation (H). This weight has traditionally been the inverse of the 
Gamow penetrability factor M. The Gamow factor is "pointlike", in that it describes 
the emission of an isolated pair of charged particles from a small and otherwise neutral 
spatial region (r < 1 fm). Because of the finite size of realistic sources [2B| and the 



multiparticle Coulomb screening which occurs when a large number of charged pions are 
simultaneously emitted ||20| , Gamow factor corrections are expected not to be appropriate 
for central Pb+Pb collisions at CERN SPS energies. 

The Coulomb correction for charged pions can be investigated experimentally through 
measurement of the correlation functions of oppositely charged pairs, where quantum me- 



chanical symmetrization effects do not contribute to the correlation function |20| . Fig. |2], 
panel (a), shows the correlation function for pairs of oppositely charged pions, obtained 
from the correlation of oppositely charged hadrons by applying a Monte-Carlo correction 
for non-pionic contaminations. They are plotted as a function of 



Q mv = ^{Ap x f + (A Py y + {Ap z f - (AE) 2 (4) 

near mid-rapidity and for a restricted K± range for central collisions of S+Ag and Pb+Pb. 
Also shown is the calculated correlation function for the point-like Gamow function (solid 
line), and the correlation function due to Coulomb attraction based upon finite source 
sizes of 4 fm (dotted line) and 6 fm (dashed line), following the approach of Baym and 
Braun-Munzinger ||26|| . It is seen that the Coulomb attraction calculated for a point-like 
source does not describe the measured data, and that the finite source size needs to be 
taken into account. 

Fig. |2|, panels (b) and (c), show the measured two-particle correlation function of neg- 
ative pions, extracted from the negative hadrons with the same procedure mentioned 
above, as a function of Q inv for Pb+Pb collisions, both uncorrected (open circles) and 
corrected for Coulomb repulsion effects (filled circles). Panel (b) uses the Gamow factor 
for a point-like source whereas panel (c) employs a parametrisation of the opposite-charge 
correlation function in panel (a). The two correction procedures generate different corre- 
lation functions, both in the peak width at low Qmv and in the asymptotic convergence 
to unity at high Qi nv . We find that the finite-size correction leads to a constant value of 
the correlation function at Qi nv > 60 MeV/c (see inset of panel (c)), while the "standard" 
Gamow correction results in a slope of the correlation function at large Qi nv and leads to 
a systematic underestimate of the source parameters for large dense pion sources^. All 



1 Note that the Gamow function for pairs of like charged particles in a very good approximation is equal 
to the inverse Gamow function for pairs of unlike charged particles; a slight difference is seen only for 
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correlation functions discussed in the remainder of this paper have been corrected for 
Coulomb repulsion using the finite-size source correction as in Fig. |2|, panel (c), obtained 
from the opposite-charge correlation function. 

4. Determination of Source Parameters from Correlation Functions 

The data set used for the two particle correlation analysis consists of 80,000 central 
Pb+Pb events. After tracking quality cuts and non-primary vertex track rejection the 
average multiplicity for a central Pb+Pb collision within 2.9 < y n < 5.5 is about 200 
negative hadrons. The resulting number of close pairs in this sample enables the analysis 
of correlation functions in the three relative momentum projections Qo, Q± and Q\\, in 
bins in K± of width 100 MeV/c and bins in pair rapidity Y n7T of width 0.5, with small 
statistical error. An independent analysis of the data of VTPC2 and the MTPCs is 
presented here, with the two data sets overlapping in the region 3.4 < y n < 4.4, where the 
results agree within ±5%. In the overlapping region the data was not averaged, but the 
low rapidity points at y n < 3.9 are from VTPC2 and the remainder are from the MTPCs. 
Preliminary results of such a differential analysis of the correlation functions have been 
presented in [^7] . 



The three dimensional fit to the expression (|2]) is performed using the maximum like- 
lihood method to fit its five parameters, with the correlation functions evaluated in the 
Local Centre of Mass System which refers to the mean rapidity in each rapidity bin. 
Fig. |3| shows a set of projections for a typical correlation function of negative hadrons 
as a function of Qj_, Q\\ and Qo, within 3.4 < < 3.9 and 0.1 < K± < 0.2 GeV/c. 
Also shown is the result of the fit with the parameters R± = 5.7 fm, Rn = 7.4 fm and 
Rq = 3.5 fm. Non Gaussian source functions have been recently proposed, e.g. 
but we find that a Gaussian source function fits the observed correlation functions at all 
rapidity and K± bins very well, as shown in Fig. |3|. 

Fig. [| shows typical fit results for the parameters R\\ , R± and Ro as a function of pair 
rapidity Y n7T within 0.1 < K± < 0.2 GeV/c and as a function of K± within the interval 
3.9 < Y wn < 4.4. At mid-rapidity (Y W7T = y cm = 2.9), R\\ and R± peak at 6-7 fm. 
The source appears to be locally isotropic. At higher rapidities the source parameters 
are reduced to about 5.5 fm. The parameter Rq, which reflects the duration of emission 
is approximately constant at Rq = At = 3 — 4 fm/c. Previous measurements of 



At in two-particle correlation studies at the SPS with Sulphur beams have yielded very 
small values of Ro, implying a "sudden" freezeout. 

We defer analysis of the correlation intensity parameter A to a forthcoming publication 
but note here that throughout this analysis, in which the pion pair purity is 55%, the fit 
result for A is about 0.45. 

Fig. |] shows the dependence of the "Yano-Koonin-Podgoretskii" rapidity Yykp as a 
function of pair rapidity Y nn for central Pb+Pb collisions. Yykp is derived from (3yk 
as 

1 1 + Pyk . 

Yykp = V y cm (5) 

z 1 — Pyk 

Qinv < 10 MeV/c. The same is true also for the model of Baym and Braun-Munzinger independent 
of the assumed particle source size. 
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Figure 2. Panel (a): Correlation function of pairs of oppositely charged pions as a 
function of Qi nv , for central collisions of S+Ag (stars) and Pb+Pb (circles), within the 
indicated acceptance. Also shown are calculated correlation functions based upon the 
Coulomb interaction for a point-like source (solid line) and finite-sized sources of radii 4 
fm (dotted line) and 6 fm (dashed line) p6| . Panel (b): Correlation function of pairs of 
negatively charged pions as a function of Qi nv , for central collisions of Pb+Pb, within the 
same acceptance as panel (a). Open circles are uncorrected data, filled circles are data 
corrected for Coulomb repulsion calculated with standard Gamow function for point-like 
source. The insert shows the same distribution with expanded vertical scale. Panel (c): 
same as panel (b), but correction now based upon finite-sized source of radius 6 fm p6| . 
The correlation functions have been derived from positive and negative hadron data by 
an appropriate Monte-Carlo correction for non-pionic contaminations. 
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Figure 3. Two-particle correlation functions for negative hadrons from central Pb+Pb 
collisions, as a function of Q±, Q\\ and Q , in the interval 3.4 < Y wn < 3.9 and 0.1 < 
K± < 0.2 GeV/c. For each panel, the correlation function is integrated from to 50 
MeV/c along the axes orthogonal to the one plotted. Solid curves are the result of a 
3-dimensional maximum likelihood fit of equation (|2]) to the data (see text). 
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Figure 4. Examples of parameters R\\ , , and Rq resulting from a fit of equation 
(|2|) to the two particle negative hadron correlation functions for central Pb+Pb collisions. 
Left panels: parameters as a function of Y^ n within 0.1 < Kj_ < 0.2 GeV/c. Filled 
circles are measured data, open circles are measured data reflected about mid-rapidity 
(Y 7Tn = 2.9). Right panels: parameters as a function of K± within 3.9 < Y n7T < 4.4. 
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where (3yk is determined by fitting expression (H) to the negative hadron two-particle 
correlation functions in five successive bins of pion pair rapidity Y n7T . This representation 
was first used by the GIBS Collaboration [^9| . Fig. || shows that the rapidity Y YKP of 



the local source is strongly correlated with the rapidity Y wn of the emitted pion pairs. 
A non-expanding source would exhibit no correlation, whereas an infinite boost-invariant 
source would exhibit a strict correlation. The data are consistent with a boost-invariant 
source, but show some deviation from the boost-invariant scenario at high rapidity, as is 
expected for a beam of finite energy (y^^ = 5.8). At Y n7r = 5 the longitudinal expansion 
velocity is found to be (3yk = 0.9. We turn next to an analysis of the transverse source 



expansion, employing the dependence of R± on K± which was predicted by Heinz [^2 



to be sensitive to the presence and magnitude of a radial velocity field. Fig. |6| shows the 
dependence of R± on K± for four different bins of Y nn , where R± is determined by the 
fit of (0) to the negative hadron two-particle correlation functions. The solid line is the 
result of a fit to the data by the function 



R± = Ro[l + cosh(Y YKP - Y m )]- 1/2 ; M T = y'm 2 + K\ (6) 

This expression was derived by Heinz |§ in a model with longitudinal and transverse 
expansion, under the assumption that the radial profile of the transverse velocity is given 
by Ptrans(r) = f3±(r/Ra), where Rq is the geometrical RMS radius of a Gaussian ap- 
proximation to the one dimensional transverse projection of the source density profile at 
freezeout, and r is the distance from the centre 0. The fit parameters are Rq and the 
ratio [3j_/T with T the freezeout temperature. The dotted lines in Fig. [] show equation 
(0) with parameters one standard deviation away from those at minimum \ 2 ■ The data 
in the three lowest rapidity bins are best described with a common set of fit parameters: 
Rq = (6.5 ±0.5) fm and f3j_/T = (3.7 ± 1.6) GeV -1 . At forward rapidity both parameters 
decrease, indicating that the pion emitting source is ellipsoidal in shape. We note here 
that the effect of long lived resonance decays on the Kx dependence is ignored in the 
model leading to @. Wiedemann and Heinz |28| have shown, however, that this effect, 



in itself produces a slope of R± with K± which leads to a relatively small apparent value 
of pl/T = 0.75 GeV" 1 if © is used to fit the R x slope. 

5. Discussion of Source Parameters 

The fit of equation ([]) to the K± dependence of R± in Fig. ^| fixes the genuinely ge- 
ometric transverse source radius Rq = 6.5 fm near mid-rapidity. In order to assess this 
result we recall the geometrical aspects of a Pb+Pb collision. We start from the 208 Pb 
Woods-Saxon density distribution of radius i?( 208 Pb) = 1.15x(208) 1 / 3 fm = 6.8 fm. In an 
ideal zero impact parameter collision the primordial participant baryons would be posi- 
tioned in a cylindrical volume of this dimension. However the NA49 trigger for central 
events selects a mean impact parameter of about 2 fm 0, resulting in a slight reduction 



2 In 1 28 1 one notices that j3 is replaced by a transverse rapidity r\ in order to be valid in the relativistic 
regime. In the domain of /3 we find in this study the differences between the two is of the order of our 
systematic error. 
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Figure 5. Yano-Koonin-Podgoretskii rapidity Yykp vs. laboratory rapidity of the 
pion pair Y T7r for central collisions of Pb+Pb. Yykp is derived from f3 Y K (equation |). 
Filled circles are data, open circles are data reflected about mid-rapidity. The dashed 
horizontal line indicates the correlation expected for a non-expanding source, the dotted 
line indicates the correlation expected for an infinite boost-invariant source. 
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Figure 6. Dependence of R± on K± in four successive rapidity intervals, for central 
Pb+Pb collisions. Solid curves show result of the fit using equation (H). Dashed curves 
show the variation of the fit function when the (3\/T parameter is varied one standard 
deviation from the best fit values. 
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of the azimuthally averaged radius, to about R par tici P ants = 6 fm. It is important, however, 
to consider next the transverse density distribution in the primordial source. Right after 
interpenetration, before significant transverse expansion has occurred, this is a longitudi- 
nal cylinder. The transverse density profile is obtained by projecting the two spherically 
uniform nucleon distributions onto a plane perpendicular to the beam (z) axis. In order to 
compare to the geometrical source parameter obtained from the correlation analysis, we 
recall that the transverse radius Rq = 6.5 fm refers to a projection of the density profile 
of the cylinder onto one transverse dimension (x) only. Applying this projection also to 
the primordial source cylinder, we obtain a nearly Gaussian density profile. Overall, the 
two projections transform the shape of two uniform spheres (where the radius is identical 
to the "surface" ) to approximately a Gaussian, with RMS radius smaller by v5 than the 
initial sphere radius. I.e. in the terminology employed for the source at freezeout, with 
RMS radius Rq = 6.5 fm as obtained above, the primordial source has an RMS radius of 
Rprim « 6/y/E fm = 2.7 fm. The source thus expands by a factor of about 2.4 between 
initial maximum energy density, and final hadronic decoupling. If it is assumed that the 
hadronic transverse expansion velocity at the RMS radius position is (3 ~ 0.5 (see below), 
the system requires a Gaussian mean expansion time of about r exp = 8 fm/c. 

The Kx dependence of the source radius parameters shown in Fig. £| exhibits a steep 
decrease in Rn and a shallower decrease of Rx- The duration-of-emission parameter Rq 
stays roughly constant with perhaps some indication of a decrease at high Kx- Heinz || 
has provided an estimate of R\\ in the framework of the same source model that leads to 
the estimate of Rx employed above (equation (^j)): 



R\\ = r 



M T 11 

— cosh (Y YKP -Y m ) -j— — — + -— 

T cosh (Y YK p - ) (A77)' 



(7) 



where r is the overall time of expansion, T is the freezeout temperature which we take 
as 120 MeV and Arj is the width of the source in space-time rapidity which we take 
to be 1.3 (see bellow). We obtain r pa 8 fm/c near mid-rapidity, decreasing slightly to 
t pa 6 fm/c at high rapidity. Note that the overall expansion time r of the system can be 
longer than the interval At during which hadrons decouple (freeze out) from the source 
(the duration of the emission); indeed we observe that r > At = Rq, where Rq is seen in 
Fig. |]to be 3-4 fm/c. Recalling the Gaussian parametrization of the correlation functions, 
the temporal emission pattern of the hadronic fireball can thus be pictured as a Gaussian 
with mean r = 8 fm/c and o = R = At = 3.5 fm/c. This is consistent with our above 
estimate of T exp , the net duration of transverse expansion. The emission of pions fades 
away at about 3a, i.e. after 8 + 3x3.5 fm/c ~18 fm/c. 



6. Transverse Expansion and Single Particle Spectra 

The Kx dependence of Rx extracted from the two-particle correlation data and the fits 
using equation (|6]) were shown in Fig. |6] for four rapidity bins. The first three rapidity 
bins are best described with a common set of parameters, Rq = (6.5 ± 0.5) fm and 
0x1 T = (3.7 ± 1.6) GeV -1 . Because the fit parameter /3j_/T is a ratio, there is an 
ambiguity in the transverse velocity fix and freezeout temperature T, ranging from fix 
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= 0.86 at T = 200 MeV to f3± = 0.61 at T = 100 MeV (see Fig. |7|). The implied 
average transverse energy of a single pion varies dramatically over this domain. We 
remove this (T,f3±) ambiguity and enforce transverse energy conservation by fitting the 
NA49 transverse mass spectra [|16|] within the same model. 
It has been shown 
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that the spectra of pions and heavier hadrons in Pb+Pb 
collisions are consistent with the assumption of a radial velocity field. We employ here a 
parametrization given by Chapman et al. [^] to fit the single particle transverse mass 
distribution for central Pb+Pb collisions. It refers to the space-time source function || as 
was employed in our above analysis of pion correlations (equation @. The single particle 
distributions in this case can be approximated by the expression 



Pi(m T ,y) oc m T RlAr]* 



* exp 

where 
1 1 



rriT 



+ 



ft 



1 + 



[Win 



3_ 

ml 



(A77): 



(A v )i 



m 



8T 
0.5y 2 



(8) 



T 2T(T + m T (3l) (Ar]) 2 + T/m T 



1+P 



2 mr 
T 



and 



(A V )1 [A V f T 



I rriT . , 

+ -^7- with irtT 



mi 



Pt 



(9) 



As already mentioned above in the context of equation (^]), the model's parameters are 
the temperature T, the transverse geometrical source radius Rq, the transverse expansion 
velocity (3±_ at the distance R and the width Ar) of the source in the space-time rapidity. 
The latter is related to the width of the rapidity distribution Ay by 



(Ay) 2 = (At?) 2 + — 
rriT 



(10) 



For negative hadrons in a realistic temperature range T from 80 to 180 MeV, Ay = 1.4 
and < rriT >= 0.45 GeV giving Arj = 1.3 ± 0.1. 

The fit results in again an ambiguous set of solutions with regard to T and /3±, but 
now with a different constraint than the two pion correlation, namely the observed mean 
transverse momentum. For negative hadrons at mid-rapidity we thus obtain "true" tem- 
peratures at freezeout ranging from T w 180 MeV at (3± = to T w 120 MeV at /3± ~ 0.7. 

The allowed regions in the plane of freezeout temperature T versus transverse surface 
velocity /3± that result both from the correlation and spectral analysis of negative hadrons 
are shown in Fig. [7[ Overlayed on the same figure is the analogous region of allowed T 
versus /3j_ obtained by NA49 from the measured transverse mass spectrum of deuterons 
|33|j which were identified by dE/dx and time-of-flight. The contribution of the background 
under the deuteron peak was determined in small bins of total and transverse momentum 
[ P4j] and has been subtracted. The bands shown correspond to a one a variation from the 
best fit values taking into account the systematic errors. For negative hadrons the choice 
of the fit interval gives a contribution to the systematic error which accounts for half of 
the width of the band, for the deuterons this is negligible. A similar analysis of (more 



preliminary) NA49 data was carried out by Kampfer |pq| . 
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Figure 7. Allowed regions of freezeout temperature vs. radial velocity for central Pb+Pb 
collisions near mid-rapidity, derived from the fit of equation @ to negative hadron two- 
particle correlation functions and from negative hadron and deuteron transverse mass 
spectra fitted by expression (§). Bands are drawn at ±cr around fitted values. 



All data employed in Fig. |^ refer to a rapidity domain about 0.6 units away from mid- 
rapidity. The negative hadron spectra and correlation data refer to 3.4 < y < 3.9, the 
deuteron spectra to 2 < y < 2.4 (kinematically equivalent relative to y cm = 2.9). The 
details of the analysis of hadronic spectra will be given in a forthcoming publication |33| . 

From Fig. [7| we see that the three independent results favour a narrowly defined overlap 
region, T = (120 ± 12) MeV and /3± = 0.55 ± 0.12. The errors given here contain 
an estimate of the systematic error of the negative hadron and deuteron spectral data 
included in the width of the bands in Fig. [?]. The remaining systematic error from the 
correlation data was estimated to be negligible for T and to be about (+0.02, -0.05) for 
/3j_. No correction has been made for the decay of long-lived resonances. However the 
deuteron data should be fairly free of such effects, and Heinz has shown [22] that the K± 
dependence of R± is only weakly affected by this at the relatively high radial velocity 
extracted here. 



18 



7. Conclusions 

We have carried out a high statistics analysis of Bose Einstein correlations of negative 
pions in central Pb+Pb collisions at 158 GeV per nucleon, which has allowed a differential 
study of the dependence of the correlation function on the pair kinematic variables, in 
small bins of rapidity and transverse momentum of the pair. The correlation between 
opposite-charge hadrons measured simultaneously in the same apparatus is used to ex- 
perimentally determine the Coulomb correction, a method that is shown to be superior 
to the correction with the standard Gamow factor. 

The representation of the correlation function derived by Heinz ^,22| which is most suit- 



able for a source with longitudinal collective motion is employed to extract the parameters 
for an expanding source. The five fitted parameters A, R±, Rn, i? an d Pyk are functions 
of the pair rapidity Y nn and transverse momentum K± . The duration of emission given 
by Rq is 3 to 4 fm/c showing that there is a finite duration of the hadronic freezeout. It is 
nevertheless small and not consistent with a long duration of a mixed (hadronic-partonic) 
phase implied by a strong first order phase transition [[/J. This statement might have to 
be modified in a scenario with significant source expansion before hadronization as was 
pointed out already by Pratt |7|]. The rapidity dependence of (3yk is found to correspond 
to a source with large longitudinal expansion, close to the boost-invariant limit. 

We then employ a model which contains longitudinal and transverse expansion || 
to extract information about the transverse dynamics of the source from the Y nn and 
K± dependence of the transverse radius R±. The results correspond to a source with 
a Gaussian radius Rq = (6.5 ± 0.5) fm, which is a factor of 2.4 larger than the initial 
transverse RMS radius of the high density zone. The transverse dynamics are indicated 
by the parameter /3j_/T = (3.7 ± 1.6) GeV -1 , which however cannot fix the transverse 
velocity and the temperature of the system independently. The use of the transverse 
momentum spectra of deuterons and negative hadrons give further constraints fixing T = 
(120 ± 12) MeV and /3± = 0.55 ± 0.12. 

With the temperature of the source fixed, the rapidity dependence of R\\ gives the 
overall expansion time r = 8 fm/c. This agrees with the expansion time needed for the 
system to grow transversely by a factor of 2.4 with a transverse velocity /3±= 0.55, showing 
the consistency of the results of the longitudinal and transverse expansion modes. 

The observation of longitudinal expansion is not unique to central nuclear collisions, and 
follows from the longitudinal dynamics of incompletely stopped partons and the decay of 
their "strings" predominantly ordered in beam direction |36j. However, no trivial expla- 
nation is at hand for the origin(s) of the large transverse collective velocity field. It may 
arise from the combined effects of the initial hadronic and partonic scattering cascades 
[ |37]j (in the extreme case from the transverse expansion of a QCD partonic state prior to 
hadronization), or from an isentropic expansion of the tightly packed hadrons following 
hadronization, that either generates a radially ordered velocity field, or enhances such an 
initial field if it exists already at hadronization |T3|J14] . 
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